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ABSTRACT 

Following the recent theoretical interpretation of GRB 090618 and GRB 101023, we here interpret 
GRB 970828 in terms of a double episode emission: the first episode, observed in the first 40 s of the 
emission, is interpreted as the proto-black-hole emission; the second episode, observed after to +50 s, as 
a canonical gamma ray burst. The transition between the two episodes marks the black hole formation. 

The characteristics of the real GRB, in the second episode, are an energy of E^ ot e = 1.60 x 10 53 erg, 
a baryon load of B = 7 x 10 -3 and a bulk Lorentz factor at transparency of T = 142.5. The clear 
analogy with GRB 090618 would require also in GRB 970828 the presence of a possible supernova. We 
also infer that the GRB exploded in an environment with a large average particle density (n) ~ 10 3 
part/cm 3 and dense clouds characterized by typical dimensions of (4 — 8) x 10 14 cm and Sn/n oc 10. 
Such an environment is in line with the observed large column density absorption, which might have 
darkened both the supernova emission and the GRB optical afterglow. 



1. INTRODUCTION 

GRB 97 0828 is one of the first GRBs with an observed 
afterglow (jPiro et al.l Il999f ) and a determined redshift 
of z=0.9578 from the identification of its host galaxy 
([Djorgovski et al.ll2Q01h . This source is still presenting 
today, after 15 years from its observations, an extremely 
rich problematic in the identification of its astrophysical 
nature. 

The X-ray afterglow was discovered by the ASCA satel- 
lite 1 .17 days after the GRB trigger (Mura kami et al.l 
Il997f ) and the data showed the presence of a fea- 
ture around ~ 5 keV, maybe associated to a radia- 
tive recombination e dge of a H- or He-like ionized iron 
([Yoshida et al.ll2001[ ). Particularly interesting was also a 
variable and large intrinsic absorption column which was 
found from the AS CA: the absence of an optical afterglow 
([Groot et al.l Fl998), the larg e intrinsic absorptio n column 
detected in the X-ray data (| Yoshida et al.ll200lt ) and the 
contemporary detection in radio- wavelengths of the GRB 
afterglow, imply a very large value for the circum-burst 
medium (CBM), this variable absorption might be an 
i ndication of a stron g inhomogeneous CBM distribution. 

iPe'er et all (|2007| ) described the presence of an evolv- 
ing thermal component in the first 40 s of the emis- 
sion of GRB 970828, as observed by the BATSE de- 
tector on board the Compton Gamma Ray Telescope 
(M eegan et all 11992( 1. This thermal emission was asso- 
ciated to the photosp heric GRB emis sion of a relativis- 
ts expanding fireball (Meszaros 2002), and from the ob- 
served properties of the blackbody spectrum it was in- 
ferred the bulk Lorentz Gamma factor of the expanding 
plasma, T « (305 ± 28). This first thermal component 
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is followed by a second emission episode, which starts ~ 
50 s after the GRB trigge r. This second ep isode was ne- 
glected in the analysis of iPe'er et al.l (2007), where they 
explicitly stated to "neglect here late-time episodes of 
engine activity that occur after ~ 25 and ~ 60 s in this 
burst" . 

We have rec ently proposed that some G RBs, namely 
GRB 09 0618 dRuffinil l20lfl Ezzo et al.ll20Tl ) and GRB 
101023 ([Penacchioni et all 120121 ), can be described by 
two sharply different emission episodes. A first episode 
that is characterized by a thermal component in the early 
phases of its emission, which originates from an expand- 
ing source, not related to the GRB emission. We have 
shown that this component is characterized by an ex- 
pansion of the blackbody emitter, whose origin is asso- 
ciated to the latest phase of the life of a massive star. 
In the late phase of its evolution, this massive star is 
well-described by a bare core that has lost its Hydrogen 
envelope and most or entirely the Helium one, during 
the previous stages of its evolution. This core is then 
observed to expand at non-relativistic velocities. At the 
end of the first episode, the collapse to a black hole oc- 
curs and the actual GRB emission happens. Since this 
first episode occurs prior to the GRB emission and to the 
formation of a black hole, we have referred to this as the 
proto-black-hole ( PBH - irpuro^ from the ancient greek 
that means anteriority in space and in time ). 

In this work we describe first the PBH component that 
we find in the first episode of GRB 970828 light curve. 
We then proceed to the description of the second episode, 
which in this work represents the authentic GRB emis- 
sion. The GRB emission is well explained in the co ntext 
of the fireshell scenario, see e.g. iRuffini et al.l ([20071 ) for 
a complete review of the model. We find for the GRB 
an e + e _ -plasma total energy of E^ t e = 1.60 x 10 53 
erg, a baryon load B — 7 x 10 -3 and a Lorentz Gamma 
factor at transparency of T = 142.5. From the numeri- 
cal simulations we confirm the presence in the CBM of 
large average particle density, of the order of (n) « 10 3 
part /cm 3 , and dense clouds with Sn/n oc 10. We show 
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Figure 1. BATSE-LAD light curve of GRB 970828 in the 25-1900 
keV energy range. 

that this can explain the "dark" nature of this GRB and 
the observed large and variable column absorption in the 
ASCA X-ray afterglow. These results are in very good 
agreement with the analysis presented in Yos hida et al.l 
(2001) and represents a further confirmation of the pres- 
ence of a PBH in the early emission of GRB 970828. 

2. DATA ANALYSIS 

GRB 970828 was discovered with the All-Sky Mon- 
itor (ASM) on board the Rossi X-Ray Timing Ex - 
plorer (RXTE) on 1997 August 28th ([Smith et al.lll997D . 
Within 3.6 hr the RXTE/PCA scanned the region of 
the sky around the error box of the ASM bur s t and 
detected a weak X-r ay source (Ma rshall et al.l 119971 : 
iFilippenko et all 11997( 1. GRB 970828 was also observed 
by the Burst and Transient Source Ex periment (BATSE ) 
and the GRB experiment on Ulysses (Smi th et al.l fl997). 
The BATSE-LAD light curve is characterized by two 
main emission phenomena, see FigJTJ the first lasts about 
40 s and is well described by two main pulses, the second 
one is more irregular, being composed by several sharp 
pulses, lasting other 40 s. The X-ray afterglow obser- 
vations started at about 0.1 days after the trigger time 
and continued up to 7-10 days. It is well described by 
a double power- law behavior, but a possible resurgence 
of soft X-ray emission, together with marginal evidence 
for the exist ence of Fe lines, wa s reported by Beppo-SAX 
and ASCA (jVietri et al.lll999f ). 

The optical observations, which started about 4 hr af- 
ter the burst, did not report any pos sible optical after - 
glow for GRB 970828 up to #=23.8 (jGroot et al.lll998ft . 
However, the observations at radio wavelengths of the 
burst position, 3.5 hrs after the initial burst, succeeded 
in identifying a source at a good significance level of 4.5 
a ([Dj orgovski et al. 2001) inside the ROSAT error circle 
(10"). The following deep searches for a possible optical 
counterpart of this radio source led to the identification 
of an interacting system of faint galaxies, successively 
recognized as the host galaxy of GRB 970828. The spec- 
troscopic observations of the brightest of this system of 
galaxies led to the identification of their redshift, being 
z=0.9578. The lack of an optical transient associated 
with the afterglow of GRB 970828 can be explained as 
due to the presence of strong absorption, due to dusty 
clouds in the burst site environment, whose presence does 
not affect the X-ray and the radio observations of the 
GRB afterglow. 

To analyze in detail this GRB, we have considered 



the observations of the BATSE-LAD detector, which ob- 
served GRB 970828 in the 25-1900 keV energy range, 
and then we have reduced the data by using the RM- 
FIT software package. For the spectral analysis we have 
considered the High Energy Resolution Burst (HERB) 
data, which consist of 128 separate high energy resolu- 
tion spectra stored during the burst emission. The light 
curve was obtained by using the Medium Energy Res- 
olution (MER) data, which consist of 4.096 16-channel 
spectra summed from triggered detectors. 

2.1. The PBH emission in the first emission episode 

In analogy to the cases of GRB 090618 (llzzo et al.l 
I2012f ) and GRB 101023 (jPenacchioni et al.ll2012l ). we an- 
alyze here the first emission episode in GRB 970828 to 
seek for a PBH signature. An integrated spectral analy- 
sis of these initial 40 s of emission is reported in Table [U 
This first time interval is best fitted by a Band spectral 
model and by a double blackbody with an extra power- 
law component. A single blackbody with power-law is 
not satisfactory, providing af = 2.10, see Fig. [2] 

A detailed time-resolved investigation of these initial 
40 s of emission has provided very interesting results: 
the analysis with a single blackbody plus a power-law 
model has revealed a double decay trend of the temper- 
ature corresponding to the two main spikes in the 
observed light curve of this first episode, see Fig. 03 We 
have then analyzed this characteristic evolution of the 
blackbody in both time intervals, corresponding each one 
to an observed decay trend of the temperature. From 
the observed flux of the blackbody component (j)BB,obs 
for each interval, we obtain the evolution of the emitter 
radius, which is given by 

__ f 0BB,obs \ D n s 

Tem ~ V oT* ) (1 + z) 2 ' 1 j 

where a is the Stefan constant, T the observed temper- 
ature, D the luminosity distance of the GRB and z the 
redshift. The evolution of the emitter radius for the com- 
plete PBH emission is shown in Fig. HI It is very inter- 
esting that the radius monotonically increases, without 
showing the double trend observed for the temperature, 
see Fig. [3l The global evolution of the emitter radius 
is well-described with a power-law function r oct 6 and a 
best fit of the data provides for the S index the value of 
0.42 ± 0.05, with an R 2 statistic value of 0.98, see Fig. 

HI 

We conclude that the double blackbody observed in the 
integrated emission of the first episode is indeed due to 
the expansion, during the PBH emission, from an initial 
radius of 7000 km to a final distances of ~ 20000 km 
from the progenitor. 

2.2. The GRB emission in the second emission episode 

The observed fluence in GRB 970828 and the redshift 
z = 0.9578 imply an isotropic energy for the total GRB 
emission Ei so = 4.2 x 10 53 erg. Following the separa- 
tion of this GRB in two different episodes we have com- 
puted the isotropic energies emitted in both episodes, 
by considering a Band model as the best fit for the ob- 
served integrated spectra: E is0y i st = 2.6 x 10 53 erg and 
Eiso,2nd = 1.6 x 10 53 erg. In what follows we explain 
the second emission episode of GRB 970828 as a single 
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Table 1 

Spectral analysis (25 keV - 1.94 MeV) of the first 40 s of emission in GRB 970828. 



Spectral model a 


P 


E (keV) 


kTi{keV) 


kT 2 (keV) 


7 


x 2 


Band -0.62 ± 0.03 


-2.14 ± 0.05 


358.0 ± 12.2 








1.21 


BB + po 






63.00 ± 0.89 




1.51 ± 0.01 


2.10 


DoubleBB + po 






108.0 ± 6.5 


40.6 ± 2.0 


1.54 ± 0.02 


0.99 
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Figure 2. Time- integrated spectral fits and sigma residual plot (25-1900 keV) of the first emission episode (0-40 s) in GRB 970828 with 
respectively a) a blackbody plus an extra power-law model; b) a double blackbody plus an extra power-law component; c) a Band model. 

canonical GRB emission in the context of the Fireshell 
scenario. 

In this mo dej_jRuffini et all l2QQlri iRuffinil [200H 
Ruffini et al.l l2Q01aL a GRB originates from an 
optically thick e + e~-plasma created in the pro- 
cess of vacuum polarization by a charged black 
hole (jChristodoulou fc Ruffinil Il97ll iDamour fc Ruffinil 
|1975[ ) . The dynamics of this expanding plasma is mainly 

described by its total energy E^ t e , the baryon load 

B = M B c 2 lE e to f and the circum burst matter (CBM) 
distribution around the burst site. The GRB light curve 
emission is characterized by a first brief emission, named 
the proper GRB or P-GRB, originating in the process of 
the transparency emission of the e + e~-plasma, followed 
by a multi-wavelength emission due to the collisions of 
the residual accelerated baryons and leptons with the 
CBM. This latter emission is assumed in a fully radiative 
regime and is called the extended afterglow. In the spher- 
ically symmetric approximation the interaction of the ac- 
celerated plasma with the CBM can be described by the 
matter density distribution ucbm around the burst site 
and the fireshell surface filling factor 1Z = A e ff/A visi 
which is the ratio between the effective emitting area 
and the total one (Ruffini "et al.ll2QQ5h . The spectral en- 
ergy distribution is well-described by a "m odified" ther- 
mal emission model (jPatricelli et al.ll2011f ). which differs 
from a classical blackbody model by the presence of a 
tail in the low-energy range. 

In this context, to simulate the second episode of GRB 
970828, which is the actual GRB emission, we need to 
identify the P-GRB signature in the early second episode 
light curve. From the identification of the P-GRB ther- 
mal signature, and the consequent determination of the 
energy emitted at transparency, we can obtain the value 
of the baryon load B assuming that the total energy of 
the e + e~ -plasma is given by the isotropic energy E iso 
observed for the second episode of GRB 970828, as it 
was done for the second episode in GRB 090618, see e.g. 
Ilzzo et all (|2012[ ). We have then started to seek for any 



Figure 3. The evolution of the temperature kT (blue bars) as 
obtained from a time-resolved spectral analysis of the first 40 s of 
emission of GRB 970828. The light curve of the first episode (red 
dots) is shown in background. 
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Figure 4. Evolution of the radius of the first episode of GRB 
970828. The line corresponds to the best fit of this dataset with a 
power-law function r oc t* 5 , with 5 = 0.42 zb 0.05. 
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thermal feature attributable to the P-GRB emission in 
the early emission of the second episode. As it is shown 
in Fig. [6j the early emission of the second episode is 
characterized by an intense spike, anticipated by a weak 
emission of 7 s. The emission observed in this time in- 
terval can be fitted by a power-law with an exponential 
cut-off model (Comptonization model in RMFIT) as well 
as with the addition of a blackbody component, see Table 
Eland Fig. El 

The observed fluence in the P-GRB emission, com- 
puted from the fit with the blackbody and the power- 
law with exponential cut-off model is S Q bs = 1.16 x 10 -6 . 
This value of the fluence is at the boundary of the region 
of minimum fluence that a burst should have in order to 
perform a spectral analysis, in a given energy range, and 
then derive the main spectral paramet ers, as the peak 
energ y or the blackbody temperature (jGhirlanda et al.l 
2008), see also Fig. [71 In this light, we can conclude that 
the signature of the P-GRB emission is represented by 
these early 7 s of emission and in the following we will 
base on this time interval for the P-GRB emission of the 
second episode of GRB 970828. 

From the observed fluence and the values obtained 
from the spectral fits, we have computed the energy emit- 
ted in the P-GRB by using the Schaefer formula of the K- 
correction (Schaefer 2007) and the redshift of the source: 
it is Epgrb = 1.46 x 10 51 erg. This value of the P-GRB 
energy is only 9 x 10 -3 the isotropic energy of the sec- 
ond episode, which we can assume as equal to the energy 
E^ ot e of the e + e _ -plasma of the GRB emission. In this 
way, with the knowledge of the temperature at trans- 
parency and of the P-GRB energy, we can derive the 
value of the baryon load from the numerical solutions 
of the fireshell equations of motion. These solutions for 
four different values of the total e+e~ -p lasma energy are 
shown in the Fig. 4 of llzzo et all (|2012[ ). We find that 
the baryon load is B = 7 x 10 -3 , which corresponds to a 
Lorentz gamma factor at transparency T = 142.5. The 
GRB emission was simulated with very good approxi- 
mation by using a density mask characterized by an ir- 
regular behavior: all the spikes correspond to spherical 
clouds with a large particle density (n) ~ 10 3 part/cm 3 , 
and with radius of the order of (4 — 8) x 10 14 cm, see Fig. 
[U Considering all the clouds found in our analysis, the 
average density of the CBM medium is (n) = 3.4 x 10 3 
particles/cm 3 . The corresponding masses of the blobs 
are of the order of 10 22 g, in agreement with the clumps 
found in GRB 090618. 

The results of the fireshell simulation are shown in Ta- 
ble [3] while the simulated light curve and spectrum are 
shown in Fig. [9] 

2.3. Limits on a possible SN-related observations and 
neo-neutron star 

The analys is of GRB 090618 (|Izmel_aL| [Soil) and 
GRB 101023 ( Penacc hioni et al.ll2012[ ) represents an au- 
thentic "Rosetta Stone" for the understanding of the 
GRB-SN phenomenon. The possible presence of a su- 
pernova emis sion, observed te n days after the burst in 
GRB 090618 (jCano et al.ll2011[ ) but not spectroscopically 
confirmed, is at the same time very intriguin g in light 
of th e still unknown GRB-SN connection (|Della Valid 
120111 ). We have transposed the data of the "bump" Rc- 



band light curve observed in the optical afterglow of GRB 
090618, associ ated to the presen ce of a possible underly- 
ing supernova (|Cano et al.H20TTl) . to the redshift of GRB 
970828. This simple operation concerns only the trans- 
formation of the flux observed, under the assumption 
that the SN has the same intrinsic luminosity: the vari- 
ation of the magnitude only depends on the square of 
the difference between the luminosity distances. More- 
over, we have also transposed the U and R— band light 
curves of SN 1998bw (?), which is the prototype of a su- 
pernova associated to a GRB. From the K— correction 
transformation formula, the U— band light curve, trans- 
posed at z=0.9578, corresponds approximately to the ob- 
served R— band light curve, so in principle we should 
consider the U = 365 nm transposed light curve as the 
actual one observed with the Rc = 647 nm optical fil- 
ter. These transposed light curves are shown in Fig. [TUJ 
We conclude that the possible Supernova emission could 
have been seen between 20 and 40 days after the GRB 
trigger, neglecting any possible intrinsic extinction. The 
optical observations were made up to 7 d ays from the 
GRB trigger, reaching a limit of R ~ 23.8 (fGroot et al.l 
1998) and subsequent deeper images after ~ 60 days 
([Djorgovski et al.ll2001[ ), so there are no observations in 
this time interval. However, the i?-band extinction value 
should be large since the observed column density from 
the X-ray observations of the GRB afterglow is large as 
well (jYoshid a et al.ll2001[ ): the computed light curve for 
the possible SN of GRB 970828 should be lowered by 
more than 1 magnitude, leading to a SN bump below 
the R = 25.2 limit, see Fig. [TUJ The presence of very 
dense clouds of matter near the burst site might be the 
cause of these absorptions. Indeed we find the presence 
of clouds in our simulation at the average distances of 
~ 10 15-16 cm from the GRB progenitor, with average 
density of (n) ~ 10 3 part/cm 3 and typical dimensions of 
(4 - 8) x 10 14 cm, see Fig. 

The possibility to observe the energy distribution from 
a GRB in a very wide energy range, thanks to the newly 
space missions dedicated, has allowed to definitely con- 
firm the presence of two separate emission episodes in 
GRBs associated to SNe. Future missions, as the pro- 
posed LOFT 0, will allow to observe the thermal decay 
from these objects down to kT = 0.5 — 1 keV. The Large 
Area Detector, designed for the LOFT mission, could be 
also able to observe in timing mode the emission from 
the newly born neutron star (NS) from Supernovae Ib/c, 
catc hing a possible chara cteristic signature of the neo- 
NS {Ne ereiros et all 12012( 1. shedding more light on the 
SN-GRB enigma. 

3. CONCLUSIONS 

A general conclusion is therefore appropriate: it is 
mandatory to proceed to a detailed spectral analysis in 
parallel to an identification of the physical and astrophys- 
ical nature of the source. This has bee n made explicitl y 
clear in the analysis of GRB 090618 (jlzzo et all 120121 ) . 
The crucial role of GRB 090618 has been to evidence the 
existence of a new family of GRBs, associated to SNe, 
composed of two distinct emission episodes. 

This has been made possible thanks to different rea- 
sons: 1) the large number of detectors, Fermi, Swift 

6 http:/ /www.isdc.unige.ch/loft/index.php/the-loft-mission 
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Figure 5. Time-integrated spectral fits and sigma residual plot (25-1900 keV) of the P-GRB emission in the second episode (0-40 s) 
in GRB 970828 with respectively a) a power-law model with an exponential cut-off; b) a blackbody plus a power-law model with an 
exponential cut-off. 

Table 2 

Spectral analysis (25 keV - 1.94 MeV) of the P-GRB emission in the second episode of GRB 970828. 



Spectral model 


7 


E cutoff (keV) 


kT(keV) 


C-STAT/dof 


Compt 


-0.68 ± 0.42 


969.8 ± 180 




75.3/11 


BB + Compt 


-1.13 ± 0.58 


894.4 ± 152 


9.23 ± 2.46 


66.7/9 
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Figure 6. Light curve of the second episode in GRB 970828. The 
dashed region represents the P-GRB emission. 



Table 3 

Final results of the simulation of GRB 970828 in the fireshell 
scenario 



Figure 7. Observer frame plane of the peak energy and bolomet- 
ric fluence (1 - 10 4 keV). The blue region represents the value of 
the observed fluence for the 7 s of emission of the second episode in 
GRB 970828, where the P-GRB emission is identified. The shaded 
orange region represent the minimum fluence that a burst detected 
by BATSE should have in order to perform the spectral analysis 
and make constraints about th e main spectral parameters. The 
figure is taken from the work of Ghirlanda et al. ( 2008). 



Parameter 



B 
To 

kT th 

Ep-GRB,th 

<n> 
5n/n 



Value 



(1.60 ± 0.03) x 10 53 erg 
(7.00 ± 0.55) x KT 3 
142.5 ± 57 
(7.4 ± 1.3) keV 
(1.46 ± 0.43) x 10 51 erg 
3.4 x 10 3 part/ 'cm 3 
10 part /cm 3 



BAT and XRT, AGILE, Suzaku-WAM, Konus-WIND 
and Coronas PHOTON RT-2, that allowed to study that 
GRB in a very wide electromagnetic energy range, from 
0.3 keV to 40 MeV; 2) the very large intensity emitted 
from GRB 090618, E iso = 2.7 x 10 53 erg; 3) the low 
redshift of the GRB, z = 0.54. Although both episodes 
radiate in the hard X-ray energy range, the first episode 
originates from a non-relativistic expansion process, hav- 
ing an expansion velocity of ~ 1000 km/s, see Fig. |4j and 
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Figure 8. Radial CBM density distribution for GRB 970828. The 
characteristic masses of each cloud are on the order of ~ 10 22 g 
and 10 15 cm in radii. The black line corresponds to the average 
value for the particle density. 
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Figure 9. Light curve with the fireshell simulation (top panel) and 
spectrum (low panel) with the fireshell simulation of the second 
episode emission in GRB 970828, in the context of the fireshell 
scenario. The red lines correspond to the fireshell simulation. 



R-band transposed light curve of SN1 998bw 
Rc-band transposed light curve of the SN-GRB 090618 
U-band transposed light curve of SN1 998bw 
R-band limit given by Groot et al. 1 999 
R magnitude of the B component described in Djorgovski et al. 2001 




Time (days from burst) 

Figure 10. The transposed data at the redshift of GRB 970828 
of the same SN data (green data) with the U— (blue line) and 
R— band (red line) light curve of SN 1998bw transposed at the 
same redshift of GRB 970828 . The purp le and cyan line repr esents 
the limit given in t he deep images by [Djorgovski et al. (2001]) and 
Groot et al. (1998) respectively. 

processes. 

We verify in this work that indeed GRB 970828 is a 
member of this family. The new understanding leads to 
a wealth of information on the evolution of the thermal 
component and of the radius of the blackbody emitter, 
given by Eq. (1), in the first emission episode, see Figs. 
[3l [U Similarly, for the second episode we obtained the 
CBM structure, see Fig. [5J the details of the simulation 
of the lightcurve and the spectrum of the real GRB emis- 
sion, see Fig. [9] We have also shown in Table [2] the final 
results of the GRB simulation, the total energy of the 
e + e~ plasma, the baryon load £?, the temperature of the 
P-GRB kT t h and the Lorentz Gamma factor at trans- 
parency T, as well the average value of the CBM density 
(^cbm) and the density ratio of the clouds Sn/n. Fi- 
nally, from the analogy with the optical bump observed 
in GRB 0906 18, see Fig. fTUl associated to a possible 
SN emission (Ca no et al.l 1201 1[ ), we have given reasons 
why a possible SN associated to GRB 970828 would not 
have been observable due to the large interstellar local 
absorption, in agreement with the large column d ensity 
observed in the ASCA X-ray data lYoshida et al.l ([2001) 
and with the large value we have inferred for the CBM 
density distribution, (ucbm) ~ 10 3 particles/cm 3 . 
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